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We investigate the effect of fc-space broadening of the interlayer pairing kernel on the critical tem- 
perature T c and the fc-dependence of the gap function in a one-dimensional version of the interlayer 
pair-tunneling model of high-T c superconductivity. We consider constant as well as fc-dependent in- 
tralayer pairing kernels. We find that the sensitivity to fc-space broadening is larger the smaller the 
width of the peak of the Fermi-level gap calculated for zero broadening. This width increases with 
the overall magnitude of the interlayer tunneling matrix element, and decreases with the bandwidth 
of the single-electron intralayer excitation spectrum. The width also increases as the Fermi level is 
moved towards regions where the excitation spectrum flattens out. We argue that our qualitative 
conclusions are valid also for a two-dimensional model. This indicates that at or close to half-filling 
in two dimensions, when the Fermi-surface gap for zero broadening attains its peaks at (±7r/a,0) 
and (0, ±7r/a) where the excitation spectrum is flat, these peaks should be fairly robust to moderate 
momentum broadening. 



I. INTRODUCTION 

The interlayer pair-tunneling (JLT) model of high-temperature supercoaductivityEI has been the focus of much 
attention since it was introducedcl and later elaborated on quantitatively.cl Within the ILT model, the pairing of 
electrons in individual Cu02-layers is considerably enhanced by the tunneling of Cooper pairs between neighbouring 
layers, giving critical temperatures which are substantially higher than those arising solely from a reasonable in-plane 
effective electron-electron attraction in a two-dimensional (2D) BCS-like theory. 

The central underlying assumption of the ILT model is that the normal state of the cuprates is a strongly correlated 
non-Fermi liquid, where sing/e-clcctron interlayer tunneling is incoherent or strongly damped, resulting in a frustrated 
c-axis kinetic energy. This kinetic energy is substantially lowered in the superconducting state, through tunneling of 
Cooper pairs between CuC"2 layers. Thus, contrary to the situation in conventional superconductors, it is the lowering 
of the kinetic energy, and not the potential energy, which drives the transition. 

Recently, there have been extensive.. discussions in the literature about experimental tests of an unconventional 
relation, predicted by the ILT model,ut3 between the c-axis penetration depth A c and the conelensation energy i5 conc j. 
The agreement in LSCO seems quite good,EfEl but experiments on Hg-120lQ and Tl-220ll§LrEij siye estimates for A c 
which are 8-20 times larger than the predicted valuesB However, note that Chakravarty et a/.llil have argued that 
this discrepancy between theory and experiment can be drastically reduced by taking more carefully into account the 
fluctuation contributions to the normal state specific heat when estimating -E C ond- 

It is not our purpose here to consider the microscopic foundations of the ILT mechanism. Instead, we will take it 
as a phenomenological starting point, and explore the effects of some modifications of the form of the pair tunneling 
term used in Ref. y. There it was argued that, in order to obtain critical temperatures of the same order of magnitude 
as found in the high-T c cuprates, it was essential that the 2D momentum of the Cooper-pair electrons was conserved 
in the tunneling process. This momentum conservation was argued to follow from the momentum conservation of the 
single-electron tunneling Hamiltonian, in the absence of inelastic scattering. Translated to real space, this momentum 
conservation means that the electron-electron attraction associated with the interlayer tunneling has an infinite range. 

A natural question to ask is then how sensitive the critical temperature T c is to a relaxation of this constraint. 
Specifically, what is the typical order of the interaction range below which T c will drop to values which are no longer 
comparable to critical temperatures in the high-T c cuprates? Moreover, several of the unusual fc-space features of the 
gap predicted within the ILT mechanism have their origin in the assumed momentum conservation. 

We will address this question phenomenologically, modelling the finite range by postulating modified functional 
forms of the pair tunneling term in which phenomenological parameters are introduced to measure the degree of 
"screening" . This will in turn lead to modifications of the original gap equations, which are then solved self-consistcntly 
to obtain the critical temperature and the superconducting gap. We expect that qualitatively correct conclusions may 
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be drawn from our modelling of the k-space broadening. Brief accounts of parts of this work have appeared in print 
elsewhere. 



II. FORMULATION OF THE PROBLEM 

For simplicity, we consider compounds with two CuQ2-layers per unit cell. The generalization to an arbitrary 
number of CuC>2-planes per unit cell is straightforward £j Below the superconducting transition temperature, we will 
assume that the quasi-particle description is approximately valid. The total Hamiltonian is taken to be the sum of 2D 
BCS Hamiltonians for the individual layers, and an interlayer pair tunneling Hamiltonian, H = i?i a ycr + Hj. When 
the zero-momentum pairing assumption is invoked, the intralayer part is given by 

n layer — ^ fc fc L fc, CT c k,a T 2-*/ k ' k k.T "M -k'.J. fc',T> ' 
fc,o-,i=l,2 fe,fc',i=l,2 

while the interlayer pair-tunneling contribution to the Hamiltonian is given by the form 

Hj = - £ T 7 (fc, k')cWc%c% lC %\ + h.c. (2) 
fc.fc' 

Here c& is the creation operator of an electron in layer i (i = 1,2) with 2D in-plane wave vector fc and spin projection 
<7, £fc is the normal state dispersion measured relative to the Fermi level, and V kt k' is the inplane contribution to the 
pairing kernel. 

An apparently pathological aspect of a particular version of (^), namely with a fc-diagonal tunneling term Tj(fc, fc') = 
Tj5k,k' p becomes evident on Fourier-transforming back to real space, where it takes the form 

J± V c (1)t c (1)t c (2) c (2) +hc (3) 

Ri,R 2 ,r 

where N is the number of lattice sites per layer, and r is the relative coordinate and Ri the center of mass coordinate 
in layer i of the two tunneling electrons . Note that there are no restrictions on \Ri — i?2| due to the zero-momentum 
pairing assumption, as in conventional superconductors. What is not conjwptional is that there is no restriction on the 
relative positions in each plane for which two electrons feel an attraction.^ Hence, TjSk,k' represents an infinite-range 
attraction, contrary to the conventional case where it is a (retarded) contact-attraction. That such a version of the 
ILT-model then gives a large value of T c is perhaps not surprising, but it is difficult to understand how such an 
effective attraction is produced. 

The fc-diagonal model must therefore be viewed as an idealization, and the issue to adress is how representative 
this limit is, if at all. The more general model given in (||) yields 

~~/\f ^2 G (l r l) C fl\ ) +r/2,T C fl\ ) -r/2 : l C fl 2 ) -r/2a C it 2 ) +r/2,T + h - C - ( 4 ) 
Rl,R2,r 

The characteristic decay-length of the function G(|r|) = J2k elkr f{k), with /(fc) defined via Tj(k, fc') = Tjf(k — fc'), 
represents the range of the effective interlayer tunneling attraction. 

By assuming a layer-independent pair amplitude, the total Hamiltonian becomes decoupled in the layer indices, and 
the gap equation is seen to be the same as in the BCS case when one makes the replacement Vk,w ~^ Vk.k' — Tj(k, fc'), 
i.e. 

A fe = ~Y^ v k,k'^k'Xk' + ^2Tj(k,k')A k , X k', (5) 
k' fc' 

where Afe is the gap function, and Xfc is the pair susceptibility, given by Xfc = tanh(/3i?fc/2)/2£ , fe, where Ef. — 
\fz\ + |Afe| 2 , (3 — 1/fcgT, fcs is Boltzmann's constant and T is the temperature. We will consider Vk,k' to be a 
separable function of fc and fc', i.e. V^k 1 = — ^5fcfffc' ; where g k belongs to the set of basis functions for irreducible 
representations of the point group of the underlying lattice, and V > is an effective two-particle scattering matrix 
element. 



2 



III. GAP EQUATION IN ENERGY SPACE 



Ref. [| studied the case Tj(k,k') — Tj6k,k', i.e. the pair tunneling matrix element is both diagonal and k- 
indepcndcnt. Using also the BCS approximation gk = Q(u>d — \ £ k\)i where u>d is an eaergy cutoff, the gap then 
depends on k only through Eu, so that the gap equation can be written in energy space asc3 

A(e) = A e(u D - \e\) + TjA(e) X (e), (6) 

where 

A = A deA(e) X (s). (7) 

The BCS coupling constant is A = VN(ep), where N(sf) is the density of states per spin at the Fermi level Ef (i-e. 
here we have made the usual approximation of neglecting the variation of V^-e inside the thin Debye shell around the 
Fermi energy) . This gap equation can be regarded as the limit u> — > of the more general equation 

A(e) = A Q(u; D - \e\) + g de' A(e') X (e'), (8) 

where the parameter lj provides a measure of the amount of fc-space broadening in the interlayer pairing kernel. 

We have solved (g) self-consistently and show in Fig. [l]the results for T c as function of w for Tj = 30 meV, ojd = 20 
meV and A = 0.1. The most important feature of this figure is the moderate reduction of T c as u> is increased from 
zero. To reduce T c by a factor 2 requires a broadening of lo ~ 40 meV. 

If we convert the energy broadening of the ILT term to a length using w = h 2 k 2 / (2M), with M equal to the electron 
mass, we obtain for the length I = l/k 

where u> is to be measured in meV. Setting uj — 40 gives an interaction range / w 9.8 A. 



IV. GAP EQUATION IN ID if-SPACE 

In this section, we will consider the gap equation (||) with a particular choice of Tj(k,k'). The main purpose of 
this paper is to establish a qualitative criterion for how robust the sharp /c-space structures of the gap, obtained for 
a fc-diagonal ILT term, are to momentum broadening. Given this limited purpose, it does make sense to simplify 
the problem by taking the fc's to be one-dimensional (ID). This simplification is purely mathematical, and of course 
does not imply anything about superconductivity with true off-diagonal long-range order in ID systems, which is 
well-known not to exist for T > JlH and prohibited by quantum fluctuations at T = 0. The final justification of our 
ID model lies in the qualitative conclusions established at the end of this section, which will be seen to apply also to 
a 2D system. 

We will consider two different functional forms for gk, both giving a fc-symmetric gap as required for singlet pairing. 



The first is the BCS approximation gk — 0(lud — \sk\), also used in Sec. HI. It is analogous to isotropic s-wave pairing 
in 2D. The second form is gk = cos(fca), which is most closely analogous to 5^2+^2 or d x 2_ y 2 pairing in 2D. The gap 
obtained for the first form does not change sign in the Brillouin zone, while the gap for the second form in general 
does. 

For simplicity, we assume a simple tight-binding dispersion form for 

Sk = —2t [cos(fca) — cos(kpa)] , (10) 

where t is the single-electron intralayer tunneling matrix element, a is the lattice constant and kF is the Fermi wave 
vector. The pair tunneling term is taken to be of the form Tj(k,k') = Tjf(k — k'), where we have chosen f(k) to 
have the particular form 

2L sin 2(^) + (Y) 
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where L is the length of the system and fco is a measure of the width of /(fc). The prefactor in ( pd| ) is chosen to ensure a 
fc-diagonal ILT term in (j^) in the limit ko — » 0. The sine function ensures that the scattering is periodic in the reciprocal 
lattice. One could construct infinitely many functions f(k) which reduce to a delta function as fco — ¥ 0, and hence 
our particular choice (|Tl| ) is inevitably somewhat arbitrary. However, since our focus here is merely on the qualitative 
aspects of momentum broadening, the detailed form of f{k) is of no concern to us; any function /(fc) which is "smeared 
out" as fco increases, would give the same qualitative results. Note that G(r = 0) = 53fc/(^) = + (k a/2) 2 , 

which means that the effective value of Tj actually decreases as fco a is increased. In this respect, the effect of 
momentum broadening is at least not underestimated in our model. 



Results and discussion 



We have calculated the critical temperature T c and the zero-temperaturegap function for various values of fco<x by 
solving (^) self-consistently in the thermodynamic limit L — > oo. In Fig. |2j we show the results for T c for Tj = 30 
meV, u>d = 20 meV, t — 25 meV, kpa — 7r/4 and LV/2-kcl = 2.5 meV. It is seen that T c is slightly more sensitive 
to fcoa for gk = cos(fca) than for g^ = 0(u)d — For g^ = cos(fca), T c is reduced by a factor 2 compared to the 

fc-diagonal result when k a/ir s» 0.25. Only 1/10 of this broadening is required for a 50% reduction of T c if one instead 
chooses Tj = 50 meV, t = 250 meV and LV/2-ira = 25 meV.EJ The reason for this increased sensitivity to broadening 
is the large increase of t. 

In Fig. U we show the gap at T — for four values of Tj and fixed fco = 0, the other parameter values being the 
same as used for Fig. ^[ In this case, the gap is given implicitly by 

= _^off^ ( } 

1-TjXk y ' 

where Ao = V^2 k gk&kXk- The maximum of the gap, and hence the critical temperature T c , is determined by Tj 
through the enhancement factor 1/(1 — TjXk), which has its maximum on the Fermi surface. However, as seen in 
Fig. ||, Tj does not affect the sign of the gap, which is determined by g^ alone. On a 2D square lattice, the analogous 
statement is that the transformation properties of the gap function under the symmetry operations of th e ffo fnt group 
of the square lattice, C± v , is given entirely in terms of the intralayer contribution to the pairing kernel,£3lI3 which is 
expandable in terms of basis functions for the irreducible representations of C± v . 

In Fig. H we show the gap at T = for four values of fcoa and fixed Tj — 30 meV. Note how the fc-space variation 
of the gap decreases with increasing fcoa. For large enough fcoa, /(fc) is essentially independent of fc, so the ILT term 
in (|J) essentially becomes a constant self-consistent shift of A/.. The main contribution to the shift comes from the 
Fermi surface region, where A& and \k are maximal. Therefore, the sign of the shift is essentially determined by 
the sign of A^ on the Fermi surface, which in turn is determined by the sign of g^ on the Fermi surface, which for 
gu = cos(fca) changes at half- filling. Thus, the qualitative form of the gap is given by A^ = Aos^ + Tj Ai, where, 
for gj~ = cos(fca), the sign of Ai is positive below half- filling and negative above half- filling. As a consequence of this 
shift, A*, eventually ceases to change sign in the Brillouin zone for g^ — cos(fca), as seen in Fig. |^. 

We now discuss the criterion for how much broadening is needed to obtain a substantial reduction of the maximum 
value of the gap, thereby smoothing out the sharp fc-space structures obtained in the fc-diagonal case. For this purpose, 
it is instructive to consider how a slightly broadened Tj(k, fc') affects the maximum value of the gap. For fc a/7r <C 1, 
Afc varies more rapidly in the Fermi surface region than Xk, because 1/(1 — TjXk) is sharply peaked at the Fermi 
surface. Thus the variation of A^x^. in the Fermi surface region is essentially determined by A&. Furthermore, the 
main contributions to Tj{kp, k^Ak'Xk' roughly come from the region \kp — k'\ ^ fco- Temporarily denoting the 
gap calculated for fco = as Afe(0), it follows that as long as fco is much smaller than the characteristic width of the 
peak of Afe(0), the broadened Tj(kp,k') essentially has the same effect as a 5-function. Under such circumstances, 
the gap is little affected by the non-fc-diagonality. A broadening of the order of the width of the peak of A/-(0) is 
therefore required for a substantial effect of the broadening to be felt. Fig. |3| shows that the width of the peak of 
Afc(0) increases with Tj. The detrimental effects on the gap of an increase of fcoa will therefore be reduced with an 
increase of Tj. On the other hand, increasing t will make the width of the peak of Afc(0) smaller, because the factor 
1/(1 — TjXk) drops more abruptly away from its peak value as one moves away from the Fermi surface when the 
overall amplitude of the variation of Sk is increased, as seen from the fact that this drop is proportional to 

fc/t = 2ta sin(fcj? a)5k. (13) 

Thus the parameters Tj and t have opposite effects on the sensitivity of the gap to broadening of Tj(fc, fc'). Note that 
one may scale the parameter t entirely out of (|^) to obtain a gap equation in terms of the dimensionless parameters 



4 



1/Pt, Tj/t, Afc/t, V/t and k^a (and uJo/t, when gt — 0(uirj — |efc|)). It should also be mentioned that the 'realistic' 
values of Tj/t are difficult to ascertain, because the model we have considered is one-dimensional, and because the 
experimentally relevant values of Tj are hard to extract. For these reasons, we can only draw qualitative conclusions 
from our model. 

Another interesting consequence of ( |l3| ) is that the width of the peak of A&(0) will increase as the Fermi level is 
moved towards the band edges where the dispersion flattens out, since then Ssk decreases. So in our ID model the 
system becomes more robust to a finite fco for a nearly empty or nearly full conduction band. 

We finally stress that our qualitative conclusions regarding the sensitivity to momentum broadening are valid also 
for a 2D model. This is because the arguments used to arrive at these conclusions depend on premises that will be 
present also in 2D: 1) in the fc-diagonal case, the gap shows sharp enhancement at the Fermi surface, 2) the width (and 
height) of the peak of this gap is increased by increasing the amplitude Tj of the interlayer tunneling matrix element, 
3) there will be parameters in the 2D single-electron intralayer dispersion analogous to t in the ID case, which control 
the bandwidth of the dispersion and therefore affect the width of the peak of the fc-diagonal gap in a manner 
similar to what occurs in the ID case. Note also that in the 2D case, the tight-binding dispersion flattens out near 
the points (±7r/a, 0) and (0,±7r/a), which lie on the Fermi surface when the band is halfjfillcd (for nearest-neighbor 
hopping only) or close to half-filling (when next-nearest neighbor hopping is included) J13 These are also the points 
where the fc-diagonal gap is at its maximum for such filling factorsJj Thus it appears that near half-filling in 2D, the 
maximum value of the gap should be fairly robust to moderate momentum broadening. 



V. CONCLUSIONS 



We have considered superconductivity within the ILT mechanism in the presence of non-fc-diagonal interlayer 
tunneling. We find that the sensitivity to momentum broadening is larger the smaller the width of the peak of the gap 
obtained for fc-diagonal tunneling. This width is increased by increasing the amplitude Tj of the interlayer tunneling 
matrix element. The width is decreased by increasing the bandwidth of the single-electron intralayer dispersion. 
Finally, the width is larger at points on the Fermi surface where the dispersion is relatively flat as compared to points 
where the dispersion is steeper.E3 Although we illustrated these features by solving a model with one-dimensional 
intralayer wavevectors, these qualitative conclusions are also valid for the more experimentally relevant case of two 
dimensions. 

Several unusual properties of the superconducting state of the cuprates are given an explanation with the ILT 
mechanism. The essential feature of the ILT mechanism is the sharp fc-space structure of the gap that arises from 
an unusual enhancement factor 1/(1 — TjXk) for a k-diagonal interlayer tunneling. Conclusions based on these sharp 
structures ought therefore to be reexamined in the presence of a slightly broadened interlayer tunneling term. This 
pertains for instance to the explanation of the anomalies in the neutron scattering peaks observed in YBCO usingr-the 
ILT mechanism. In this case, non-trivial Fermi surface kinematics almost unique to the mechanism are essen.tial.Ea 
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FIG. 1. The critical temperature T c of the gap equation (|^) as a function of the energy broadening ui. 

FIG. 2. The critical temperature T c as a function of k^aj-K. 

FIG. 3. The T = gap function plotted for four values of Tj with fco = 0. The two cases g k — Q(u>d — \s-k\) an d g k = cos(fca) 
are shown in the upper and lower panel, respectively. As Tj is increased, the maximum value of the gap, occuring on the Fermi 
surface, increases and the variation of the gap with k is enhanced. Note how the sign of the gap is always determined by g k . 

FIG. 4. The T = gap function plotted for four values of koa/n with Tj = 30 meV. The two cases g k — @(lud — £fc|) and 
g k — cos(fca) are shown in the upper and lower panel, respectively. As fco is increased, the maximum value of the gap and the 
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